Diffusive Radiation in One-dimensional Langmuir Turbulence 
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We calculate spectra of radiation produced by a relativistic particle in the presence of one- 
dimensional Langmuir turbulence which might be generated by a streaming instability in the plasma, 
in particular, in the shock front or at the shock-shock interactions. The shape of the radiation spectra 
is shown to depend sensitively on the angle between the particle velocity and electric field direction. 
The radiation spectrum in the case of exactly transverse particle motion is degenerate and similar to 
that of spatially uniform Langmuir oscillations. In case of oblique propagation, the spectrum is more 
complex, it consists of a number of power-law regions and may contain a distinct high-frequency 
spectral peak. The emission process considered is relevant to various laboratory plasma settings 
and for astrophysical objects as gamma-ray bursts and collimated jets. 
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Various kinds of two-stream instability including lon- 
gitudinal, transverse, and oblique modes [3L0,Slj| are 
common for many astrophysical objects [EH, B] and lab- 
oratory plasma settings [8|, |9[ . These instabilities are effi- 
cient to produce a high level of magnetic and/or electric 
turbulence in the plasma. We note that the electric tur- 
bulence is of primary importance for many applications, 
in particular , it is capable of charged particle acceler- 
ation 0, [3, [llj in contrast to purely magnetic turbu- 
lence, which is only efficient for angular scattering of the 
charged particle while cannot change the particle energy 
directly. Remarkably, recent theoretical and numerical 
studies [H, [13, [I3| demonstrate that the energy density 
going to the electric (Langmuir) turbulence can be very 
large, exceeding, for example, the energy density of the 
initial regular magnetic field in the plasma. 

Relativistic electrons propagating in the plasma with 
developed Langmuir turbulence will random walk due to 
random Lorenz forces produced by these turbulent elec- 
tric fields. Apparently, during this random diffusive mo- 
tion the electron will generate electromagnetic emission. 
We will refer to this radiative process as Diffusive Ra- 
diation in Langmuir waves (DRL) to emphasize the key 
role of the diffusive random walk of the particle in the 
stochastic electric fields. 

Curiously, the theory of DRL has not been developed 
in sufficient detail yet, although a few particular issues 
related to this radiative process (called also Electrostatic 
Bremsstrahlung) have been considered [3, 3 12, 3 H, 
[20j . However, the detailed treatment of the DRL spectral 
shape in various regimes is currently unavailable. This 
letter attempts to remedy the situation by calculating 
the DRL spectrum within the perturbation theory and 
determines the region of applicability of the perturbative 
treatment. 

Electromagnetic emission produced by a charged par- 
ticle can be calculated within the perturbation theory 
when the particle moves almost rectilinearly with almost 



constant velocity. Apparently, the non-zero acceleration 
of the particle in the external field should be taken into 
account to obtain a non-zero radiation intensity. This 
perturbative treatment is widely used because of its sim- 
plicity. Frequently, one calculates first the particle ac- 
celeration w(t) due to a given field along the rectilin- 
ear trajectory and then uses this expression obtained for 
w(t) to find the radiation spectrum. In case of a ran- 
dom external field, however, when w(t) is also a random 
function of time t it is more convenient to express the 
radiation intensity via spatial and temporal spectrum of 
the external electric and/or magnetic field directly 21 1. 
Accordingly, the spectral and angular distribution of the 
emission produced by a single particle in a plasma with 
random field has the form: 
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where 7* = ^7 ^ + ^ j , Q, M, and 7 are the 

charge, mass, and Lorenz-factor of the emitting parti- 
cle, c is the speed of light, Upc is the plasma frequency, V 
is the volume of the emission source, Fq^^qj^ is the tem- 
poral and spatial Fourier component of the Lorenz force 
transverse to the emitting particle velocity. 
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9 is the emission angle relative to the particle velocity 
vector V, co is the frequency of the emitted wave. Contri- 
bution Wq^ (marked with the superscript ±) is provided 
by a component of the particle acceleration transverse to 
the particle velocity. In case of the electric E (in contrast 
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to magnetic) field, there is also a component of the ac- 
celeration along the particle velocity. The corresponding 
contribution has the form 
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which is typically small by a factor 7"^ compared with 
the transverse contribution. Nevertheless, there exist 
special cases when the transverse contribution is zero or 
very small and the parallel contribution comes to play. In 
particular, for the considered here one-dimensional tur- 
bulent electric field the parallel contribution will domi- 
nate for a particle moving along the field direction. 

Let us start with the case when the particle moves 
at a large angle to the Langmuir turbulence direction 
"i? 3> 7~^, so the standard transverse contribution dom- 
inates. Following the derivation given in 21|, but with 
the Lorcnz force F = QE specified by electric E in place 
of magnetic B field, it is easy to find 
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where Kai3{qo,<l) = Cai3K{qo,q), T is the total time of 
emission; CafS describes the longitudinal nature of the 
Langmuir waves, i.e., Cap = qaqp/q'^, while K{qQ,q) is 
the temporal and spatial spectrum of the Langmuir tur- 
bulence. 

The developed approach allows for arbitrary 
anisotropy of the turbulence, although we have to 
specify the shape of the turbulence spectrum to promote 
further the calculation of the DRL spectrum. Analytical 
and numerical studies of the two-stream instabilities 
suggest that the Langmuir turbulence produced is 

frequently highly anisotropic [11, S, i, El, [3, [3 , which 

is confirmed also by available in situ observations, e.g., 
performed in the Earth magnetosphere |22j. Although 
this anisotropy can be reduced at later stages of the 
nonlinear turbulence evolution due to randomization 
of the wave vector directions, here we assume that 
the Langmuir turbulence is highly anisotropic, namely, 
one-dimensional. As an example, we can suppose that 
all the wave vectors are directed along the shock normal 
n, therefore. Cap = naUp, while the spectrum K{qQ,c\) 
can be approximated by a power-law over g|| above 
certain critical value fco: 
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Here, the presence of the second 5-function is related to 
the assumption that the electric turbulence is composed 
of Langmuir waves all of which oscillate in time with 
the same frequency Wpe; the normalization constant is 
set up by the condition J K{qo,q)dqodci — (El), where 
is the mean square of the electric field in the Lang- 
muir turbulence. 

Now, substituting (|4]) with ([5]) into general expression 
([1]), taking the integrals over dqo, dq\\, and dq± with the 
use of three available i5-functions and dividing by the to- 
tal (infinite) time of emission T we find 

_ a,k^-^ (El) Q'' sin^ d /u^ " 
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where 1} is the angle between the particle velocity v 
and vector n. This expression diverges formally when 
cost? — > 0. Accordingly, for a particle moving trans- 
versely to Langmuir turbulence direction we need to re- 
calculate the integrals in Eq. ([T|) taking into account that 
S{uj' — qo + qv) 6{uj' — qo) for qv = 0, which yields 
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in full agreement with the results of obtained for spa- 
tially uniform Langmuir oscillations. This is not a ran- 
dom coincidence. Indeed, since the ID Langmuir waves 
experience spatial variations along only one direction (of 
vector n), thus, the particle moving transversely to this 
direction "feels" spatially uniform field pattern like that 
considered in [l7| . 

There is another special geometry when intensity Q is 
insufficient to describe the DRL spectrum: it is the case 
of particle motion along vector n. Indeed, the "parallel" 
contribution becomes important for ^ 7~^- Substi- 
tuting ^ into ([3]) and taking the integrals similarly to 
derivation of (El) we find: 
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Apparently, spectra ^ and ([8]) look rather differently 
compared with the spectrum ([7]) produced by a particle 
moving transversely to vector n. Given, in particular, 
that the frequency w' is directly linked with the emission 



angle 6, the (5-function 5{ijj' 



e) permits emitting a 



single frequency only in each direction. By comparison. 
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no (5-function enters ([6]) and ([8]), thus, a continuum spec- 
trum rather than distinct frequencies is emitted along any 
direction. Clearly, there remains a distinct contribution 
to the emission intensity when lo' ^ Wpe- However, the 
range of the parameter space where this resonant con- 
dition holds is relatively narrow, so the "non-resonant" 
contribution from the remaining part of the parameter 
space where uj' ^ Wpe can easily dominate the resonant 
contribution. To see this explicitly, consider the radia- 
tion intensity into the full solid angle by integration of 
^ over (Kl = sinOdOdip w 2TTd{uj' /uj) that yields 



M2cV|cos'(9| 



1/27.2 




-u/2 



(9) 



Let us analyze essential properties of the DRL spec- 
trum on the basis of the asymptotic evaluation. At low 
frequencies to ^ Wpo7^, we can discard uj' in ^ every- 
where in the braces except narrow region of parameters 
when oj' k, Wpc- This means that for lo ^ Wpc7^ the inte- 
gral is composed of two contributions. The first of them, 
a non-resonant one, arises from integration over the re- 
gion, where uj' <^ Wpc. Here, the emission is beamed 
within the characteristic emission angle of ~ along 
the particle velocity. The integral converges rapidly, and 
so it may be taken along the infinite region, which pro- 
duces a flat radiation spectrum, lui oc o;*^, or I^^ oc 
at lower frequencies, oj < Wpc7. However, as far as lo' 
approaches Wpo, a resonant contribution comes into play. 
Now, in a narrow vicinity of cjpe, we can adopt 
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which results in a single- wave-like contribution, oauj^. 
The full spectrum at w < a;pe7^, therefore, is just a sum 
of these two contributions. At high frequencies, uj :s> 
<^pe7^i the term uj' dominates in the braces, so other 
terms can be discarded. Thus, a power-law tail, I^^ oc 
uj~'^ arises in this spectral range. 

Although this integral cannot be taken analytically in 
a general case, it is easy to obtain corresponding spectra 
numerically. Figure [T] displays a numerically integrated 
example of the DRL spectra for various angles between 
V and n for highly relativistic particle with 7 = 10^. The 
black/solid curve shows the DRL spectrum arising as the 
particle moves strictly perpendicular to vector n, which 
coincide with the spectrum arising in the case of spatially 
uniform Langmuir oscillations [17|. The spectrum con- 
sists of a rising region I (x uj^ at uj < 2u;po7^ ^ind drops 
abruptly to zero at > 2ujpcj'^. Remarkably, there are 



prominent differences between this spectrum and those 
generated for oblique particle propagation even though 
the spectra can be similar to each other in the immedi- 
ate vicinity of the spectral peak. 

To find the applicability region of the perturbation the- 
ory applied above, we should estimate the characteristic 
deflection angle of the emitting electron on the emis- 
sion coherence length Ic = 2c7^/tj, where the elemen- 
tary emission pattern is formed [2l|. Consider a sim- 
ple source model consisting of uncorrelated cells with the 
size Iq = l-KcjujQ^ each of which contains coherent Lang- 
muir oscillations with the plasma frequency Wpo. Inside 
each cell the electron velocity can change by the angle 
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6*0 ^ ujst/ {ujpel) if t^o ~ Wpe, whcrc ujst = Q (Ej) /Mc 
(in the other case, ujq > Wpe, the results of [2l| apply). 
Then, after traversing N = Ic/lo cells, the mean square 
of the deflection angle is 9'^ — 9qN ~ uj'^^ujq / (ujuj^^) . The 
perturbation theory is only applicable if this diffusive de- 
flection angle is smaller than the relativistic beaming an- 
gle, 7~^, i.e., it is always valid at sufficiently high frequen- 
cies w > = Wg4Wo7^/Wpo- Note, that the bounding fre- 
quency increases with ujq, while Diffusive Synchrotron 
Radiation (DSR) in random magnetic field displays the 
opposite trend. The perturbation theory will be applica- 
ble to the entire DRL spectrum if the condition 6^ < 7"^ 
holds for the frequency Wpe7 PH . where the coherence 
length of the emission has a maximum. This happens for 
the particles whose Lorenz-factors obey the inequality 



7 < ^po/(^sVo)- 
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Let's compare the DRL spectra calculated numerically 
(Figure [I]) for the case of a broad turbulence distribu- 
tion over spatial scales with the DSR spectra in case of 
stochastic magnetic fields 2l| . Apparently, there is a re- 
markable difference between these two emission mecha- 
nisms, especially in the case of the long- wave turbulence, 
wolcosi^l = fcoc|cosz?| ^ Wpc. First, we note that the 
perturbation theory of DRL has a broader applicability 
region, in particular, it applies for higher energy electrons 



than the perturbative version of the DSR theory [21 1 
since the criterion (fTT|) is softer than the corresponding 
criterion in 21| . This happens because in the presence of 
the Langmuir turbulence the rapid temporal oscillations 
of the electric field direction substantially compensate an- 
gular defiections of the particle, so the average trajectory 
is much more similar to a straight line than for the case 
of the random magnetic field with the same ujq and ujgt- 
Then, a distinct spectral peak at uj = 2a;pe7^ is formed 
with the linear decrease of the spectrum with frequency, 
which is not present in case of the DSR. At lower frequen- 
cies, however, this falling part of the spectrum gives way 
to a flat spectrum, which is entirely missing within the 
one- wave approach 17| and absent for exactly transverse 
particle motion. Position of the corresponding turning 
point depends on the ojq cos t^/wpc ratio in such a way 
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index a up to 0. The DRL spectral asymptote (x u^, 
which appears just below the spectral peak at 2wpo7^, 
fits well to those spectra. Remarkably, the flat lower- 
frequency asymptote, I^j oc u;°, can account for the phe- 
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nomenon of the X-ray excess 2^ 24 [ and prompt optical 
flashes accompanying some GRBs. 

In addition, this mechanism along with the DSR in 
random magnetic fields 25| can be relevant to the UV-X- 
ray flattenings observed in some extragalactic jets. For 
example [i^, X-ray observations of the jet in 3C 273 
look inconsistent with the standard synchrotron of DSR 
models. Remarkably, the entire UV-to- X-ray spectrum of 
3C 273 might be produced by DRL, which can be much 
flatter than usual DSR in the range wbc7^ oj ^ i^poT^- 



FIG. 1: DRL spectra produced by a particle with 7 = 10 
in a plasma with developed Id Langmuir turbulence for var- 
ious particle propagation directions: cosi9 = 0, solid/black 
curve; cost) — 10~^, dashed/red curve; cos-d = 0.5, dash- 
dotted/green curve; cosiJ = 1, dotted/blue curve. Param- 
eters are given m the Figure. The "parallel" contribution 
(blue/dotted curve) is very small (10~^^) for the highly rela- 
tivistic particle, although it becomes competing for moderately 
relativistic particles. 



that for ujQ I cos d \ 



> 



the flat spectral region entirely 



dominates the range from cijpo7 to a;pc7 . It is worth em- 
phasizing that the deviations of the DRL spectrum from 
the single- wave spectrum (read, from the transverse case, 
cos ?9 = 0) is prominent even for oblique propagation an- 
gles only slightly different from 7r/2, e.g., cosi? = 10"'^ 
as in Figure [TJ Therefore, the presence of a broad turbu- 
lence spectrum considered here in detail results in impor- 
tant qualitative change of the emission mechanism, which 
cannot generally be reduced to a simplified treatment re- 
lying on the single-wave approximation with some rms 
value of the Langmuir electric field. 

Modern computer simulations of shock wave interac- 
tions, especially in the relativistic case, suggest that the 
energy density of the excited Langmuir turbulence can be 
very large [l3|, e.g., far in excess the energy of the initial 
regular magnetic field. In particular, at the shock wave 
front the electric field can be as strong as the correspond- 
ing wave-breaking limit, i.e., ujst ~ ^pe [IJ]- In this case 
the random walk of relativistic electrons in the stochastic 
electric field can give rise to powerful contribution in the 
nonthermal emission of an astrophysical object, entirely 
dominating full radiation spectrum or some broad part 
of it. 

Although any detailed application of the considered 
emission process is beyond the scope of this letter, we 
mention that the DRL is a promising mechanism for the 
gamma-ray bursts and extragalactic jets. In particular, 
some of the prompt gamma-ray burst emission displays 
rather hard low-energy spectra with the photon spectral 
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